BEL ROBINSON ENERGY AND CONSTANT MEAN 
CURVATURE FOLIATIONS 



LARS ANDERSSON 

Abstract. An energy estimate is proved for the Bel-Robinson energy 
along a constant mean curvature foliation in a spatially compact vacuum 
spacetime, assuming an L°° bound on the second fundamental form, and 
a bound on a spacetime version of Bel-Robinson energy. 



1. Introduction 

Let (M,g) be a 3+1 dimensional C°° maximal globally hyperbolic vac- 
uum (MGHV) space-time, which is spatially compact, i.e. M has compact 
Cauchy surfaces. One of the main conjectures (the CMC conjecture, see 
for background) concerning spatially compact MGHV spacetimes states 
that if there is a constant mean curvature (CMC) Cauchy surface Mq, in 
such a spacetime M, then there is a foliation in M of CMC Cauchy surfaces 
with mean curvatures taking on all geometrically allowed values. Specifi- 
cally, in case the Cauchy surface Mq is of Yamabe type —1 or 0, then the 
mean curvatures take all values in (— oo, 0), or (0, oo), depending on the sign 
of the mean curvature of Mq, while in case Mq is of Yamabe type +1, the 
mean curvatures take on all values in (— oo, oo). 

The only progress towards proving the CMC conjecture so far has been 
made under conditions of symmetry, cf. , or curvature bounds [3] , [H] . 

One approach to the CMC conjecture is to view it as a statement about 
the global existence problem for the Einstein vacuum field equations 

R a p = 0, (EFE) 

in the CMC time gauge. It is known that in the CMC gauge with zero shift, 
the (EFE) are non-strictly hyperbolic [8. while in other gauges such as wave 
coordinates, the (EFE) form a system of quasi-linear wave equations for the 
metric g. In this context, it has been conjectured that the Cauchy problem 
for the (EFE) is well-posed for data in H 2 x H 1 (the H 2 conjecture, see 

HDD- 

From this point of view it is interesting to consider continuation principles 
for the (EFE), in CMC gauge. In this note we will use a scaling argument 
to prove an energy estimate for CMC foliations. The energy we consider is 
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a version of the Bel-Robinson energy. For a spatial hypersurface M in M, 
the energy expression we consider is defined by 

Q(M)= [ (\E\ 2 + \B\ 2 )p g , 

where E, B are the electric and magnetic parts of the Weyl tensor (defined 
w.r.t. the timelike normal T of M. Roughly speaking, Q bounds Cauchy 
data (g, K) on M in H 2 x H . Here, g is the induced metric on M and K is 
the second fundamental form of K. Therefore, if the H 2 conjecture is true, 
apriori bounds for the Bel-Robinson energy can be expected to be relevant 
to the global existence problem for the (EFE). 

Let H = tiK denote the mean curvature, and assume the CMC gauge 
condition H = t. Define the spacetime Bel-Robinson energy of a CMC 
foliation Ti = {M t ,t £ 1} by 

Q(JFi)=fdtf N(\E\ 2 + \B\ 2 )^ g , (1.1) 

Jl JM t 

where N is the lapse function. 

We are now ready to state our main result 

Theorem 1.1. Let (M,g) be a MGHV space-time, and let I = (t—,t+) with 
— oo < t- < t+ < 0, be such that there is a CMC foliation Ti in (M,g). Let 
t G /. 

Suppose that limsup t _^ Q{t) = oo. Then at least one of the following 
holds: 

(1) limsup 4/4+ Q(^ ,t)) = oo, 

(2) limsup 4/4+ "^ff^ 00 = oo. 

The time reversed statement with t + replaced by i_ also holds. 

Remark 1.1. Let (M, 7) be a compact hyperbolic 3~manifold with sectional 
curvature —1. Then the metric 7 = —dp 2 + p 2 7 on M = (0, 00) x M 
is fiat. It follows from the work of Andersson and Moncrief |7] that for 
small perturbations of (M,j), there is a global CMC foliation J~[t ,t) ^ n the 
expanding direction, and for this foliation, the Bel-Robinson energy decays 
a>s Q(t) = 0(H 2 (t)), which implies that the space-time Bel-Robinson energy 
Q(J~) is bounded in this case. It is interesting to consider the behavior of 
Q(F[t ,t)) when to decreases. 

The proof of Theorem II. II is based on a scaling argument, which we now 
sketch. The statement of the theorem is symmetric in time, but here we 
consider only the future time direction, the argument in the reverse direction 
is similar. Suppose for a contradiction there is a constant A < 00 so that 
Q(F[t ,u)) < A > \\ K \\loo/H 2 < A for all t G [t ,Q, and that 

limsup Q(t) = 00 . 

t/u 
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An energy estimate shows that rh(t)Q(t) < C, where is an L l,p harmonic 
radius, for some fixed p, 3 < p < 6, and hence if is bounded from below 
there is nothing to prove. Suppose for a contradiction that r/, -> as 
t /* t*. The combination r^Q is scale invariant, and hence by rescaling g to 
g' = r^ 2 g, we get a sequence of metrics g' with Q bounded. Q bounds g' in 
L 2,2 and hence we may pick out a subsequence of (g',N'), which converges 
weakly to a solution (goo, ^oo) of the static vacuum Einstein equations, 

AN = 0, (1.2a) 

V 2 iV = iVRic. (1.2b) 

It follows from our assumptions that the limit g^ is complete, and the 
limiting Noo is bounded from above and below. Then by |2j, goo must be 
flat, with infinite harmonic radius, which contradicts r' h = 1, by the weak 
continuity of on L 2 ' 2 . We conclude that in fact is bounded away from 
zero, and hence that Q does not blow up, which proves the theorem. 

2. Preliminaries 

For a space-like hypersurface M in M we denote its timelike normal T 
and induced metric and second fundamental form (g,K). We assume all 
fields are C°° unless otherwise stated. Let lower case greek indices run 
over 0, . . . , 3 while lower case latin indices run over 1, . . . , 3. We work in an 
adapted frame e a , with eo = c\. Our convention for K is K a b = — \CTQ a b, so 
that if the mean curvature H = tiK is negative, T points in the expanding 
direction. We will sometimes use an index T to denote contraction with T, 
for example ut = u a T a . 

In a nonflat spatially compact, globally hyperbolic, vacuum spacetime, the 
maximum principle implies uniqueness of constant mean curvature (CMC) 
Cauchy surfaces. In particular, each x € M is contained in at most one 
CMC Cauchy surface, and for each t 6 1, there is at most one Mt with 
mean curvature t. 

Let / C R be an interval. A CMC foliation T\ in M is a foliation Ti = 
{M t ,t G 1} such that for each t £ I, M t is a C°° CMC Cauchy surface with 
mean curvature t. When convenient we will write g(t),K(t) for the data 
induced on Mj. Introducing coordinates x a with x° = t, the lapse and shift 
N, X of the foliation are defined by d% = NT + X. We may without loss of 
generality assume X = 0. 

We call Ti a maximal CMC foliation in M if there is no interval I 1 con- 
taining I as a strict subset with a CMC foliation Tj>. Given a foliation in 
M, we write Mjr for the support of T. 

Assume that M contains a compact, constant mean curvature (CMC) 
Cauchy surface Mq with mean curvature H° < 0. By standard results there 
is then an interval I = C M, H° € /, such that there is a CMC 

foliation J-j, and by uniqueness, M^o = Mq. Hence if Ti is a maximal 
CMC foliation, then I is open. 
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2.1. The Bel— Robinson energy. Let be the Weyl tensor of (M,g) 
and let *W denote its (left) dual (in vacuum *W = W*). The Bel-Robinson 
tensor Q of (M, g) is defined by 

QaM& = W*w w A V + * w ^ w ps- (2-1) 

Then Q is totally symmetric and trace-less, and in vacuum, Q has vanishing 
divergence. 

Let E a p = W a T/3T, B a/ 3 = *W a T(3T be the electric and magnetic parts of 
the Weyl tensor. Then E, B are symmetric, t-tangent (i.e. E a T = B a T = 0) 
and trace invariant, g ab E ab = g ab B a b = 0. 

In vacuum, we have 

E a b = Ric a 6 + HK ab — K ac K c b , 
B ab = -cur\K ab , 
where for a symmetric tensor in dimension 3, 

curL4 o6 = ^(e a st V t A sb + e b st V t A sa ). 

Recall that for symmetric traceless tensors in dimension 3, the Hodge system 
A I—*- (div.A, curL4) is elliptic. 

The following identities, see relate Q to E and B, 

Qtttt = E ab E ab + B ab B ab = \E\ 2 + \B\ 2 , (2.2a) 
QaTTT = 2(E A B) a , (2.2b) 

QabTT = -{E x E) ab -(Bx B) ab + ^(l^f 2 + \B\ 2 ) gij , (2.2c) 
where by definition, for symmetric tensors A, B in dimension 3, 
(A A B) a = e a bc A b d B dc , 

(A x B) ab = e a c \ ef A ce B df + ^{A ■ B)g ab - ^(tvA)(trB)g ab . 

From equation H2.2a() it follows that Qtttt ^ with equality if and only 
if W = 0. Let J- be a foliation in M. The Bel-Robinson energy Q(t) of 
Mt € J- w.r.t. the time-like normal T, is defined by 

Q(t) = Q(M t ) = [ Qtttt fJ-g = [ (\E\ 2 + \B\ 2 )^ g . 

JM t JM t 

An application of the Gauss law gives in vacuum, 

d t Q(t) = "3 I NQ^TT^Vg, 

where ir a /3 = V ' a Tp. A computation shows that the only nonzero components 
of ir a/3 are -K ab = -K ab , ir Ta = iV _1 V a iV. Thus 

NQ apT T^ = -NQ abTT K ab + Q aTT TV a N. (2.3) 
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3. Proof of Theorem 11.11 

We will assume that the complement of points ITl2*l of Theorem 11.11 holds . 
and prove from this that Q(t) does not blow up. Assume for a contradiction 
there is a constant A > 1 so that for t E [to,t*), 

GCW)) < a- ^ff < A > 

and that limsup t y tif Q(t) = oc. 

We let L s ' p denote the IP Sobolev spaces and write H s for L S)2 . We will 
sometimes use sub indices x or t,x to distinguish function spaces defined 
w.r.t. space or space-time. 

For a foliation J-j, we may without loss of generality assume that T = 
N~ 1 dt, where A > is the lapse function of the foliation. Then g is of the 
form 

g = -N 2 dt 2 + g ab dx a dx b . 
The lapse function satisfies 

-AN+ \K\ 2 N = 1, (3.2) 
which using the maximum principle implies the estimate 

l/ll-Klll- < N < 3/H 2 . (3.3) 
Let T be a foliation in M. From ()2.3j) we get 

\d t Q\ < Ci(||VJV||l» + ||iV|k-||^IU-)Q, 

with C\ a universal constant. By assumption, \K\ 2 / H 2 < A. Let K = 
K - (H/3)g be the traceless part of K. Using \K\ 2 = H 2 /3 + \K\ 2 > H 2 /3, 
we get 

\d t Q\<C(\\VN\\ Laa + A) Q . (3.4) 

Ln | 

We may assume without loss of generality that ||ViV||x,oo > A/\H\, since 
otherwise there would be nothing to prove. Therefore we may absorb A/\H\ 
in the constant in (|3.4|) to get 

\d t Q\<C\\VN\\ L ooQ. (3.5) 

3.1. The blowup. Choose once and for all a fixed p satisfying 

3<p<6. (3.6) 

On the Riemannian manifold (Mt,gt), let r^(x) denote the L 1,p harmonic 
radius of (M t ,gt) at x € M t ; thus r^(x) is the radius of the largest geodesic 
ball about x on which there is a harmonic coordinate chart in which the 
metric coefficients g a b satisfy 

rh(x)" 3/p ||5 ab -5 afe || L p( B 4 rh(x ) ) )+r/ l (x) (p_3)/p ||5g afe || LP (^( rh ( x))) < C, (3.7) 

where C is a fixed constant (say C = 1), cf. P3I1]. By Sobolev embedding, 
in B x (rh(x)), the C@ norm of g a {, is controlled, for (3 = 1 — |. The presence 
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of the factors of r^(x) in (|3.7j) means the estimate (|3.7[1 is scale invariant. It 
follows from this that rh{x) scales ELS cL distance. 

It is well known that the Laplacian in such a local harmonic coordinate 
chart on B x {rh(x)) has the form 

An = g ab d a d b u. 

Thus, within B x (rh(x)), A is given in these local coordinates as a non- 
divergence form elliptic operator, with uniform C@ control on the coefficients 
g ab , and uniform bounds on the ellipticity constants. 

We have the following standard (interior) L p elliptic estimate for this 
Laplace operator, c.f. [HJ Thm. 9.11]. Let B = B x (rh(x)) and B' = 
B x {\r h (x)). Then 

\\N\\ L 2, Pm < C(r h (x),p)[\\AN\\ LP{B) + \\N\\ LP{B) }. (3.8) 

We drop the dependence on p, since p is fixed. We need to make explicit the 
dependence of the constant C on r^(x). This is done by a standard scaling 
argument. Thus, assume (by rescaling if necessary), that r/j(x) = 1. Then 
(j3.8j) becomes 

1 1-^1 \l 2 <p(B') < C[\ \AN\\ LP ( B ) + \\N\\ LP ( B )]. 
By Sobolev embedding, since p > 3 is fixed, and B' = B(^), we have 
l|ViV|| LO o( B /) < c- \\N\\ L 2, P ( BI) , 

so that 

llwiu<» (B0 < a[||Aiv|| LP(B) + h^Hlp^)]. 

and in particular, 

||ViV|| i0 o (B) < C7 [||AJV|| L co (B) + ||iV|| L oo ( B)], (3.9) 

where C Q is an absolute constant, (i.e. independent of N, given control on 
A from definition of = 1). Now we put in scale factors to make (13.911 
scale invariant and write ()3.9|) as 

r h (x)\\VN\\ Loa{B/) < C K(x) 2 ||AiV|| LO o (B) + ||JV|| L oo (B )]. (3.10) 

Note that the function N is itself scale invariant. Each term in H3.1U|) is 
invariant under scaling, and thus (|3.1()j) holds in any scale. Therefore, it 
holds in the metric g(t). 
Let 

Th = r h (t) = inf inf r h [x). 

s<t xeM s 

From the lapse equation, 

AN = N\K\ 2 - 1. 
Using dEHIJ) and \K\ 2 = H 2 /3 + \K\ 2 we find 
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Thus, we have 



r h \\VN\\ Loa{B/) <3C (r 2 h A + ^\ . (3.11) 



In particular, this gives the estimate 

||VA^)|| L ~ <C(A,Q/r h (t). (3.12) 

Integrating ()3.5j) . (recall we absorbed the term A/\H\ in (|3.4j) into the 
constant), gives 

Q(h) < Q{t )+C I' ds\\VN(s)\\ L o°Q(s). 
Jt 

We may without loss of generality assume the last term is bigger than 1, so 
we may absorb Q(io) into C. 

Multiplying both sides by r^, and using ()3.11|) we have 

2/ 1 



r h Q(h) <C^ ds + J CM- 

We may without loss of generality assume < l/|i?|, since otherwise there 
would be nothing to prove, and therefore we can absorb the term r^A into 
the constant. Then we have 

r h Q(h)<C J*' dsj^Qis). 
The inequality ()3.3|) implies 

N > SA^H- 2 , 

which by the definition of the spacetime Bel-Robinson energy Q[F\ to ) t 1 )), 
see (|1.1|) . gives 

r h (h)Q(h) < CQ(F [t0)tl) ). (3.13) 
By assumption, limsup t y>^ Q(t) = oo, which by the assumed bound on 
2(-^[*o,ti)) im P lies 

lim rh(t) = 0. 

t/t, 

We will show that this contradicts (|3.1|) . 

Suppose then that there is an increasing sequence of times ij, U — > as 
z — ► co, so that rj = r/j(tj) satisfy lim^oo r j = (recall that by construction 
Th{t) is decreasing). 

Now we have from (|3.13|) and our assumptions, 

nQ{U) < C. (3.14) 

Now we introduce the blowup scale. Let g[ = r~ 2 g. We will denote the 
scaled versions of g, K by g^K-. We scale the coordinates as = r^t, 
x[ = r^ 1 x', so that the coordinate components of g[ are scale invariant. 
Then \K[\ < Ar^, while the lapse N does not scale, A^'(x-) = N(x). After 
translating the time coordinate as t[ = r~ l (t — ti), we focus our attention 
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on the time interval t\ G [—1,0]. We further translate the space coordinate 
so that the center of the coordinate system (0,0) is the point where the 
harmonic radius achieves its minimum value. 

Since r^Q is scale invariant, we have r^Q = Q' , and hence the inequality 

Q!{t'i) < C (3.15) 

holds. This means in view of the definition of the Bel-Robinson energy that 
at the blowup scale, Ric^ is bounded in L 2 . By construction r' h > 1 and by 
it follows from the Ricci bound that g[ is bounded in L^' . Similarly the 
Hodge system relating K to B leads to K[ bounded in L\' 2 C . 

The Einstein vacuum equation is scale invariant, and therefore holds at 
the blowup scale. We will argue in the next section, that the above bounds 
on g'ijK'i allow us to pick out a weakly convergent subsequence of (g^N^) 
with limit goo, Noq solving the static vacuum Einstein equation, cf. equation 
(|1.2|) below, with g^ complete. 

3.2. Weak convergence. Let g[ be the sequence of rescaled spacetime met- 
rics. We consider rescaled time t[ in the interval [—1,0]. By construction, 
the L 1,p harmonic radius satisfies r^(0) = 1, and r'^t) > 1 for t G [—1,0]. 
Equation (|3.3|) implies that the rescaled lapse is bounded from above and 
below, 

' <Ni<t. (3.16) 



t 2 _A ~ 1 ~ t\ 

By (I3~TB1 . we have Q!(t) < C, for t G [—1,0], and hence we have {g'^K'^ 
bounded in L°°([-l,0];Lj£ x L\£). It follows that there is a subsequence 
which converges weak-* to a limit [goo-, K^) G L°°([— 1, 0]; L 2 ^ c x L^), with 
corresponding spacetime metric goo. By passing to a further subsequence 
if necessary, which we still denote using the index i, we may assume that 
9i(0) s>oo(0) weakly in L 2 £. 

Let us consider the properties of this limit. First note that \K^\ < Ar.; — > 
as i — > oo, and hence = 0. The relation d t ig[ = —IN^K^ holds in the 
limit and since = 0, we conclude that is time independent, so that 
the limiting spacetime metric goo is static. The lapse equation now implies 
that the limiting lapse function satisfies 

Aoo^oo = 0, (3.17a) 

where Aoo is the Laplace operator defined w.r.t. g^. 

The rescaled spacetime metrics g\ are solutions of the Einstein vacuum 
equation and the evolution equation for K, 

d t K = -V 2 N + iV(Ric + HK -2K : K), 
where K : K a ^ = K ac K c b , holds weakly in the limit. In view of the fact that 

2 2 

g i— > Ric is weakly continuous on L^ c and K^o = 0, we get the equation 

= -V^iVoo + iVooRkw (3.17b) 
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By construction, is complete, and hence in view of (j3.17j) we have a 
complete solution of the static Einstein equations with Nqo > 0. It follows 
by [21 Theorem 3.2] that is flat and iVoo is constant. In particular, 
r /i[9oo](0) = oo. Now, since is by definition the L 1,p harmonic radius, 

2 2 

3 < p < 6, the map g — > is weakly continuous on L,' and hence by 
construction r/ l [5 oo ](0) = 1. This is a contradiction, and it follows that in 
fact liminfj^oo > 0, which by the BR energy estimate (|3.14j) implies that 
Q{t) does not blow up. This completes the proof of Theorem ll.il 
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